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In this article, one of the robust systems of soft computing 
namely adaptive neuro-fuzzy inference system (ANFIS) is 
used to estimate the punching shear capacity of the concrete 
column-slab connections reinforced with RP. For _ this 
purpose, a collection of experimental tests was used to train 
and test the ANFIS model. Five parameters including the 
area section of the column, Young’s modulus of the FRP 
bars, the effective flexural depth of the slab, RP 
reinforcement ratio and also the compressive strength of 
concrete are used as inputs of the neuro-fuzzy system to 
estimate the considered output. The whole structure of the 
ANFIS also presented in mathematical steps. The obtained 
results of the created model of this paper indicated that the 
proposed ANFIS structure with a suitable accuracy could be 
used as a predictive model to determine the punching shear 
capacity of the considered elements. Also, the formulated 
model of the ANFIS in this paper can easily apply for codes 
and other researches. 


1. Introduction 


One of the possible failure modes in reinforced concrete slabs is punching shear [1], and in order 
to improve the shear behavior of these types of elements, the use of fiber reinforced polymer 
(FRP) in concrete slabs is expanding. FRP is a material that can use as reinforcement materials or 
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for external strengthening of existing slabs [2]. They are brittle with a higher strength than steel 
bars and can improve the punching shear capacities of concrete slabs [3]. Theodorakopoulos and 
Swamy [4] investigated the effect of FRP on punching behavior in concrete slab-column 
connections and presented a method to design this type of elements. Zhang and Zhu [5] 
presented a verified method to finite element analysis of FRP concrete slabs, and for this 
purpose, they suggested a new element with considering nonlinearity in material and geometric 
properties. Polak and Lawler [6] studied the retrofit methods of concrete slab-column joints and 
proposed a procedure based FRP which was improved the ability of these connections in a deal 
with punching shear. Some researchers investigated the punching shear of slab-column joints 
reinforced with FRP based on experimental studies [3,7—10] and in some cases, they also 
presented models to estimate the punching capacity [7,10] or to control it [9]. Nguyen and 
Rovnak [7] proposed a fracture mechanics based on the obtained results of experimental tests of 
slabs reinforced with glass FRP to determine the punching shear capacity. They considered 
several parameters such as depth and size effect to specify the capacity. Pre-stressed FRP plates 
are considered as a suitable element in strengthening slab-column connections and the ability and 
positive effects of them on the punching shear capacity has been shown by researchers [8]. The 
effects of flexural reinforcement and space of shear reinforcement on experimental column-slab 
connections reinforced with GFRP and steel bars studied by ElGendy [9]. He was found that 
well-anchored shear studs can control the shear cracks and also the shear failure mode in these 
slabs. The influence of other parameters such as compressive strength of concrete on FRP-RC 
slab-column connection was also investigated by the researcher [10], and it was concluded that 
with increasing this parameter, the capacity of the slab is increased. A state of the art review 
regarding the punching shear capacity in slab-column joints strengthening with FRP materials 
has been presented by Saleh et al. [11]. 


The conventional methods used to determine the performance and behavior of the elements are 
very time-consuming especially in large and complex problems [12]. For this reason, in recent 
years, alternative methods such as soft computing approaches (such as neural networks, Fuzzy 
systems) has been widely used especially in civil engineering [13—26]. They are powerful tools, 
especially for capacity prediction goals. The application of artificial neural networks to 
determine the punching shear in concrete slabs reinforced with FRP bars has been investigated 
by researchers [25,27]. Despite the precision of these models, due to the complexity of the 
structure of the neural networks, the direct use of them for predictive purposes is very complex 
and time-consuming. For this reason, other soft computing methods, such as neuro-fuzzy models, 
are more useful. Such models have the benefits of neural network and fuzzy models simultaneity 
and are widely used in engineering problems because of their high accuracy and the possibility to 
provide a mathematical structure based on their solution. In this paper, the authors considered 
one of the powerful neuro-fuzzy methods namely adaptive neuro-fuzzy inference system 
(ANFIS) to estimate the punching shear of the slab-column connections reinforced with FRP 
bars. The presented model is determined based on the experimental tests and the results with 
high accuracy indicated that the ANFIS model could be used for estimating punching shear. To 
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increase the workability of the proposed model of this paper, the authors presented their model 
based on mathematical forms, and for this purpose, the equations and other requirements are also 
provided in details. 


2. Database 


Neuro-fuzzy systems are soft computing methods that need a database as an initial knowledge 
about the system behavior. These types of models are trained based on a series of data from the 
previous solutions of the problem and be able to estimate the output for an arbitrary input vector. 
In this paper, the author used a collection of 82 experimental datasets which are published in the 
literature [7,27,28]. This database was divided into two groups including 64 and 18 datasets for 
training and testing the neuro-fuzzy model, respectively. Details of the considered datasets are 
shown in Table 1. 


Table 1 
A summary on the considered database. 
Section ‘ ae Concrete Young’s FRP Punching 
Statistic aor flexural sombressive: znosulus of reinforcement — capacity 
coun depth of sucuew nS ARE ratio in percent (KN) 
(mm2) slab Gram) (MPa) slab (GPa) 
Minimum 50.27 55.00 26.00 28.40 0.18 61.00 
Maximum 2025.00 284.00 118.00 147.60 3.78 1600.00 
g Median 506.25 131.00 38.60 48.20 0.73 274.00 
= Average 762.95 138.02 42.13 70.63 0.88 455.97 
& St.Dev 575.74 58.95 14.00 35.44 0.66 351.01 
Mode 400.00 142.00 36.30 48.20 0.73 237.00 
Range 1974.73 229.00 92.00 119.20 3.60 1539.00 
Minimum 50.27 61.00 29.00 40.70 0.28 78.00 
Maximum 2025.00 284.00 75.80 113.00 3.78 1492.00 
g Median 762.50 131.00 38.60 48.15 1.12 393.00 
A Average 792.12 146.44 40.41 55.85 1.16 497.17 
& St.Dev 577.37 66.14 9.99 19.55 0.79 376.90 
Mode 900.00 131.00 38.60 42.00 1.21 N/A 
Range 1974.73 223.00 46.80 72.30 3.50 1414.00 


The authors considered five parameters as independent variables to determine the punching shear 
in concrete slab-column connections reinforced with FRP (Fig.1). These parameters are area 
section of the column, Young’s modulus of the FRP bars, the effective flexural depth of the slab, 
FRP reinforcement ratio and also the compressive strength of concrete. The considered variables 
(inputs and output) are presented in Table 2. The normal values are calculated based on the 
minimum and the range of the parameters, which can be seen in Table 1, to limit the intervals 
between 0.1 to 0.9. The use of normal value instead of real value can improve the performance of 
the final ANFIS model. 


H. Naderpour, M .Mirrashid/ Journal of Soft Computing in Civil Engineering 3-1 (2019) 16-26 


Table 2 
Definition of the considered parameters. 
Notation _ Description Normal Value 
X, 50.27 
XxX Section area of column (mm2) x, = 0.8 ——— +0.1 
1974.73 
; XA; 35 
De The effective flexural depth of slab (mm) x, = errr a +0.1 
. X;, 26 
X, Concrete compressive strength (MPa) 3 = 0.8 09 +0.1 
X, 284 
X, Young’s modulus of the FRP (GPa) , = 0.8 a6 ; 0.1 
X;, 0.18 
X; FRP reinforcement ratio in percent x, =0.8— 36 +0.1 
Punching capacity of the FRP-reinforced V 6l 
Vv ai =0.8 +0.1 
concrete flat slab (kN) 1539 
= po * 
x 8 = © 8 
= 00 i 6 $ 2 500 . a 
Ah feet er Ete 
oo moa : 
= - *63 ¢ e = ° 
S soo e ad = 500 ®—% A : 
ers tee art ee 
cog 
: 0 : rie $ 
A beh cee 


FRP reinforcement ratio in percent 


Fig. 1. Distribution of the parameters versus the punching shear. 


3. Proposed ANFIS model 


19 


In this paper, ANFIS as a neuro-fuzzy inference system, which was introduced by Jang [29], is 


considered to determine the punching shear capacities in concrete slab-column connections 
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reinforced with FRP bars. Each input variables in the proposed ANFIS structure have four 
Gaussian membership functions (MF). The selected algorithm to generate ANFIS here was 
subtractive clustering. In this method, datasets are assigned to groups called clusters to determine 
the existing patterns in the system. To estimate the unknown parameters of the ANFIS 
(membership function parameters and also the rules and the coefficients of the linear output 
functions of each rule), the model was trained based on 64 datasets and then was evaluated by 18 
separated datasets as the testing data. Details of the ANFIS is presented and formulated in the 
next section. The membership functions of the inputs are illustrated in Fig. 2. The final structure 
of the ANFIS with four rules is also shown in Fig. 3. 


Membership Functions of Input 1 ' Membership Functions of Input 2 


Membership Functions of Input 3 1 Membership Functions of Input 4 


Membership Functions of Input 5 


Fig. 2. Membership functions of the input parameters. 
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Fig. 3. The final structure of the ANFIS model. 


As previously mentioned, the presented neuro-fuzzy model of the ANFIS in this paper used four 
Gaussian membership functions (MF) for each of its input parameters. A Gaussian function can 
be shown by Eq.1 (x is the input, c is the mean and o is the variance for x). Details of the 
membership functions are presented in Table 3. Also, there were four linear functions (Eq. 2) for 
the proposed ANFIS with the parameters of Table 4. 


(xc) 


MF(x: ,c)=e 2° (1) 
Cj = A,X, +A,X, + A,X, + A,X, +asX5 + ay (2) 
Table 3 
Details of the membership functions. 
Input 1 Input 2 Input 3 Input 4 Input 5 
o c o c o c o c o fs 


MFI 0.1530 0.6754 0.0906 0.4774 0.0575 0.2843 0.0842 0.2171 0.0753 0.3060 
MF2 0.1319 0.6177 0.1105 0.5123 0.0682 0.3146 0.1737 0.7239 0.0475 0.1936 
MF3 0.1038 0.2824 0.0897 0.3234 0.0476 0.1927 0.0944 0.2159 0.0841 0.2849 
MF4 0.1103 0.1965 0.1029 0.2527 0.0455 0.2109 0.1433 0.5824 0.0696 0.1987 


Table 4 
Details of the linear functions. 
a, a, a, a, as Ay 
Cl 0.2323 1.1320 0.3688 -0.3967 0.0072 -0.2640 
C2 -0.1913 1.0020 0.4249 0.7834 1.1320 -0.7735 
C3 0.3199 0.1540 0.0009 0.2108 0.0798 -0.0188 


C4 -0.0184 0.7607 0.5423 -0.0377 0.2791 -0.1369 
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4. Mathematical steps to predict the punching shear 


In this section, the proposed ANFIS model is expressed in mathematical steps to provide a more 
efficient and simple operation of the model. The steps are as follow: 


Step 1. Determine the normal value of each input variable based on the equations of Table 2. 


Step 2. Determine the weight of each four rules based on Eqs. 3: 


(x, 0.6754) \ ( (x, 0.4774) \ ( (x3 0.2843)" \ ( (x4 0.2171) \ ( (xs 0.3060)" 


e 0.0468 } lf 0.0164 } ie 0.0066 } he 0.0142 } & 0.0113 ] 


(x, 0.6177" \ ( (x> 0.5123)" \ ( (x3 0.3146)" \ ( (x4 0.7239)" \ ( (x5 0.1936)? | 


W, = ( e 0.0348 e 0.0244 e 0.0093 e 0.0603 e 0.0085 
t \ \\ \\ \\ J 
(3) 
( _c 02824" \ (Gy 03234 \( _@y 0.1927 \( _(y_ 0.2159)? ) (_Grs_ 0.2849)? ) 
VWe\e wos e 0.0161 e 0.0045 e 0.0178 e 0.0141 
a \\ J\ \\ \\ J 
( _G 0.1965") (Gy 0.2527 \( _y 0.21097) (_ (x 0.5824) (_Gs_0.1987) ) 
W.=|e 0028 e 00212 e 0.0041 e 0.0411 e 0.0097 
=| J\ J\ }\ J\ J 
Step 3. Determine the output linear functions of each rule by Eqs. 4: 
C, = 0.2323x, +1.1320x, +0.3688x, 0.3967x, +0.0072x, 0.2640 
C,= 0.1913x, +1.0020x, +0.4249 x, + 0.7834x, +1.1320x. 0.7735 
2 1 2 3 4 5 (4) 


C, =0.3199x, +0.1540x, +0.0009x, +0.2108x, +0.0798x, 0.0188 
C,= 0.0184x, + 0.7607x, +0.5423x, 0.0377x, +0.2791x, 0.1369 


Step 4. Determine the output of the ANFIS by Eq. 5. It is worth noting that the output of the 
ANFIS is based on the normal value of the target, which is between 0.1 to 0.9. Therefore, it 
needs to convert to its corresponding real value by Eq. 6. 


WC, +W,C,+ WiC, +W,C,) 2 iC 
W,+W,+W,+W, Je 

1539Y 153.9 

+ 

0.8 


o1s[¥= 


V(kN) = 61 (6) 


H. Naderpour, M .Mirrashid/ Journal of Soft Computing in Civil Engineering 3-1 (2019) 16-26 23 


5. Results of the proposed model 


In this section, the obtained results of the proposed ANFIS are presented. The considered 
laboratory datasets contain 82 datasets, which were divided into two groups of 64 and 18 datasets 
for training and testing phases, respectively. The reason for using test data was evaluating the 
final proposed ANFIS model. The results of the training phase for all 82 datasets is presented in 
Fig. 4. It is clear from the figure that the model with suitable accuracy trained well. The amount 
of root mean square error (RMSE) and also error mean were 0.025404 and -1.86e-08, which is 
indicated that the ANFIS determined the output of 82 datasets with low error. This issue is 
illustrated in the figure for each of dataset. 


, Train Data 
Targets 
Outputs 
20.8 
= 
© 0.6 
20.4 
30 
— 
nm 
= 0.2 
0 
0 10 20 30 40 50 60 70 
Number of Data 
ei MSE = 0.00064535, RMSE = 0.025404 ¢ Error Mean = -1.8642e-08, Error StD = 0.025605 
0.05 6 
£ g 
£ 0 34 
- 3 
a = 
co 
-0.05 2 
-0.1 0 | l l 
0 10 20 30 40 50 60 70 0.1 -0.05 0 0.05 0.1 
Number of Data Errors 


Fig. 4. Results of the training datasets. 


After training the ANFIS, it should test by several new datasets. For this purpose, the authors 
used 18 datasets. The results of the proposed ANFIS structure for the testing phase is shown in 
Fig. 5 and it is indicated that the model can estimate the punching shear capacities of the test data 
with high accuracy. RMSE and error mean for these datasets were 0.04345 and 0.016371, 
respectively. For better comparison, the regression plots of the training and testing phases of the 
proposed ANFIS model is also presented in Fig. 6. It can be seen from the figure that the 
correlation coefficients for both phases were more than 0.98. 
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Test Data 


Targets 
—— Outputs 


= 


Targets and Outputs 


ie 


0 2 4 6 8 10 12 14 16 18 
Number of Data 


015 MSE = 0.0018879, RMSE = 0.04345 . Error Mean = 0.016371, Error StD = 0.041415 


Errors 
Frequency 


0 5 10 15 20 O15 4.1 0.05 0 0.05 0.1 0.15 
Number of Data Errors 


Fig. 5. Results of the testing datasets. 


Train Data: R=0.9901 Test Data: R=0.98022 


Output 


0.2 0.4 0.6 0.8 02 #03 O04 O5 06 0.7 08 
Target Target 


Fig. 6. Regression plots for the proposed ANFIS. 
6. Conclusion 


Determination of capacity of structural elements in civil engineering is an open issue that has 
widely considered in research and studies. In recent years, soft computing methods are used for 
predicting goals instead of traditional approaches such as regressions. They are simple, flexible 
and powerful and their abilities are proved in several studies. In this paper, one of this method, 
which is called ANFIS is considered to estimate the capacity of concrete column-slab 
connections reinforced with FRP. A collection of 82 experimental datasets that are published in 
literature was used to train and also test the final structure of the ANFIS. The results of the 
proposed model showed that ANFIS could be able to present the target with high accuracy. The 
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presented model considered five inputs as the independent variables and used four Gaussian 
membership function in its structure. For simplicity, the authors also presented their model in a 
mathematical form by some equations. The output results of the equations are compared with the 


experimental results, and it was concluded that the predicted values are very close to the 
experimental values. The amounts of RMSE, R2 and also error mean for both phases of training 
and testing confirm the performance of the proposed ANFIS for the considered goal. 
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